r 


AD-A071  048 


BINARY  SYSTEMS  INC  SAlIDA  CO  F/6  20/4 

FURTHER  RESULTS  ON  AN  ARRAY  DESI6N  TO  MEASURE  THE  WAVENUMBER-FR— ETC (U) 
FEB  79  L W BROOKS  N00014-72-C-0318 


UNCLASSIFIED 


CBR-79-01 


NL 


MICROCOPY  RESOLUTION  US!  CHARI 

NAIIONAl  Ht • K | At i in  -I ANI'AKi's  A 


ADA071048 


\ 


^^i7ary 


>- 

Q- 

O 

CJ> 


©LEVEL 


f 


303  E Street,  P.O.  Box  1105,  Solida,  CO  81201  (303)  539-3060 

Report  NoA  CBR-79-^l| 


FURTHER  J1ESULTS 
ON  AN  ARRAY  DESIGN 


TO  MEASURE  THE  WAVENUMBER- 
FREQUENCY  ^SPECTRUM^BENEATH 
AN  AXISYMMETRIC  ^TURBULENT 

BOUNDARY  JLAYER  . j 


1 


Final 

'7°. 


BY 


LOWELL  W. /BROOKS 

2-  /a 


5"  FebeseagR  36979/' 


D D C 

r?ror?nnnr? 

JUL  1 1 MW 

ILlbEU  UHd 


@M}~~ 


Performed  Under  ONR  Contract 
NOipjijjril  4-7  2-C- 03 18_) 

vr 


FOR 


OFFICE  OF  NAVAL  RESEARCH,  CODE  222 
800  N.  QUINCY  STREET 
ARLINGTON,  VIRGINIA  22217 

ATTENTION:  Mr.  H.  Fitzpatrick 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  release; 
Distribution  Unlimited 


^// 


1-1 


1.0  INTRODUCTION 

This  report  is  an  extension  and  correction  of 
previous  results  [1] . Here  the  requirements  on  an  array 
designed  to  measure  the  wavenumber-frequency  spectrum  of 
the  wall  pressure  beneath  a turbulent  boundary  layer  will 
be  discussed.  Of  specific  interest  will  be  the  axisymmetric 
component  of  the  pressure  field  on  a axially  towed  cylinder. 

The  basic  configuration  consists  of  an  array  of  flush- 
mounted  ring  transducers  where  the  axis  of  the  array  is 
coincident  with  that  of  the  cylinder  as  sketched  in  Figure  1.1. 
The  major  design  parameters  are  the  number  of  elements,  Ng, 
their  separation,  d,  and  width,  w.  Additionally,  the 
channel- to-channel  tolerances  on  amplitude,  phase  and 
spacing  errors  will  be  considered. 

Figure  1.2  illustrates  the  basic  problem.  There, 
two  models  which  have  been  proposed  to  describe  PQ(k,(jj), 
the  wavenumber-frequency  spectrum  of  the  axisymmetric 
component  of  wall  pressure,  are  compared  ([2],  [3]).  The 
comparison  has  been  computed  for  a frequency  of  ^ = 100  sec  * 
and  a 3 inch  diameter  cylinder  towed  at  15  kts.  The 
normalization  for  the  spectra  are  consistent  with  that  used 
by  Chase  [2].  Appendix  A gives  the  definition  of  the 
wavenumber-f requency  spectrum  which  is  used  here  and  details 
explicit  formulas  for  the  Chase  and  Gardner  models  which 
were  used  to  compute  Figure  1.2.  Symbols  which  are  not 
defined  in  the  text  can  be  found  in  the  symbol  table  at  the 
end  of  this  report. 

For  numerical  orientation.  Figure  1.2  shows  some 
wavenumber  regions  of  particular  interest.  The  acoustic 
region  corresponds  to  wavenumbers  bounded  by  |k|  < u)/c 
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Figure  1.2.  Comparison  of  Gardner  and  Chase  Models  of 
the  Wavenumber -Frequency  Spectrum  for  the  Axisymmetric 
Component  of  Pressure  under  a Turbulent  Boundary  Layer 
on  an  Axially  Towed  Cylinder. 
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where  c is  the  acoustic  propagation  velocity.  Since  the 
goal  of  a towed  array  is  to  detect  acoustic  signals,  energy 
at  these  wavenumbers  must  be  passed  with  no  attenuation. 

Thus,  the  TBL  spectral  level  here  places  an  ultimate  lower 
bound  on  the  towed  array  self  noise. 

In  the  wavenumber  reqion  for  which  |k|  = io/co, 
where  cp  is  the  breathing  (or  bulae)  wave  speed,  the  hose 
transfer  function  of  a towed  array  reaches  a maximum.  For 
a small  group  of  hydrophones  located  near  the  axis  of  a 
towed  array,  the  hose  transfer  function  governs  the  direct 
TBL  contributed  pressure  level.  Thus,  spectral  levels  in 
this  region  can  be  used  to  estimate  self  noise  for  low 
scattering  towed  array  modules. 

The  convective  region  near  k = ut/U^.,  where  Uc  is 
the  convection  velocity,  contains  most  of  the  TBL  enerqy. 

At  hose  impedance  discontinuities  such  as  endcaps,  internal 
spacers,  hydrophones  or  telemetry  cans,  energy  from  this 
region  can  be  scattered  into  hose  resonance  wavenumbers 
and  significantly  affect  towed  array  self  noise. 

The  two  models  shown  on  Figure  1.2  are  significantly 

different  at  all  wavenumbers.  The  Chase  model  indicates  a 

low  wavenumber  content  approximately  22  dB  less  than  that 

of  the  Gardner  model.  Conversely,  the  Gardner  model 

indicates  some  15  dB  less  energy  at  the  convective  peak. 

The  two  models  cross  over  at  an  intermediate  wavenumber 
k -1 

( tc—  =■  2.2  ft  in  this  case).  Although  the  numerical 

£.  IT 

values  are  different,  similar  differences  between  the  models 
occur  at  different  tow  speeds,  array  diameters,  and 
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* 

frequencies . 

Since  there  are  significant  differences  between 
the  models,  it  is  reasonable  to  ask  whether  existing 
experimental  data  can  be  used  to  indicate  which  model  is 
more  nearly  correct.  Although  a detailed  demonstration 
will  not  be  given,  it  appears  that  the  answer  is  negative. 

In  fact,  each  model  appears  to  fit  a selected  but  different 
subset  of  the  available  data. 

In  particular,  the  numerical  values  for  the  free 
parameters  in  the  Chase  model  were  obtained  by  fitting  the 
laboratory  experiments  of  Willmarth  and  his  associates  [4], 
[5].  In  contrast,  the  parameters  in  the  Gardner  model  were 
evaluated  using  a combination  of  the  Willmarth  data  and 
data  taken  at  sea  by  Markowitz  [6] . 

A major  difference  between  the  laboratory  and  the 
sea  tests  appears  in  the  frequency  power  spectrum  of  the 
pressure  fluctuations  measured  by  a small,  flush-mounted 
transducer.  Figure  1.3  compares  the  Chase  or  Gardner 
models  with  this  spectrum  as  measured  by  Willmarth  et.al. 
or  Markowitz. 

The  Chase  model  fits  the  Willmarth  data  well  for 
wa/U  < 6.  For  larger  values  of  wa/U  Chase  attributes  the 
differences  between  the  data  and  model  to  area  averaging 
by  the  pressure  transducers.  Chase  has  developed  alternate 
models  to  fit  the  Markowitz  data  which  will  not  be  discussed 
here  as  Chase  apparently  does  not  recommend  those  models. 


Both  models  (as  used  here)  can  be  written  in  the  dimensionless 
form  — P = p v*a2f(^p,  kn) ; thus  the  important  parameters 
for  model°comparisons  are  the  Strouhal  number,  wa/u  and 
normalized  wavenumber,  ka.  (See  Appendix  A) 


L . __  J 
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Area  averaging  should  not  be  a factor  in  the  Markowitz  data. 

The  Gardner  model  was  chosen  to  fit  the  Markowitz 
data  as  shown  in  Figure  1.3.  Note  that  for  wa/U  = 3.0 
(corresponding  to  Figure  1.2),  the  difference  between  models 
on  Figure  1.3  is  about  10  dB.  This  accounts  for  much  of 
the  difference  between  the  convective  peaks  seen  in  Figure  1.2. 

It  is  not  clear  (for  either  model)  that  the  same 
parameter  set  will  suffice  to  describe  both  laboratory 
and  sea  tests.  In  particular,  Willmarth  et.al.  [4]  have 
noted  that  the  symmetry  of  the  boundary  layer  is  quite 
sensitive  to  alignment  and  straightness  of  the  cylinder. 

Those  experiments  indicated  that  the  cylinder  should  be 
straight  to  within  a small  fraction  of  its  diameter  before 
axially  symmetric  flow  is  obtained.  This  criterion  may 
not  be  satisfied  by  flexible  cylinders  towed  at  sea. 

It  should  be  emphasized  that  neither  model  is 
based  on  direct  experimental  measurement  of  the  low- 
wavenumber  portion  of  the  TBL  wall  pressure  spectrum. 

The  Gardner  model  is  entirely  empirical  with  parameters 
derived  from  measurements  made  by  small  transducers.  The 
Chase  model  is  semi-empirical  in  that  the  shape  of  the 
spectrum  at  low  wavenumber  has  been  derived  from  theory, 
but  again,  the  free  parameters  have  been  derived  from 
small  transducer  measurements. 

Some  verification  of  the  models  at  low  wavenumber 
is  available  from  self  noise  measurements  made  on  low 
scattering  towed  array  modules.  In  this  case,  the  hose 
transfer  function  serves  to  reject  high  wavenumber  energy. 

If  no  scattering  occurs  near  the  hydrophone,  the  pressure 


Figure  1.3.  Comparison  of  the  Chase  and  Gardner  Models  for  the  Point  Pressure  with 
the  Experiment  Measurements  of  Willmarth  and  Markowitz. 
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measured  by  an  internal  hydrophone  is  dominated  by  the 
spectral  power  near  the  hose  resonance  wavenumbers. 

Unfortunately,  such  internal  measurements  do  not 
provide  a clear  cut  test  of  either  model.  First,  the 
models  tend  to  cross  over  near  typical  hose  resonance 
regions.  Second,  relatively  minor  modifications  in  the 
model  parameters  can  cause  significant  differences  in 
estimated  spectral  levels  (e.g.  a 20%  change  in  the  assumed 
ratio  of  v*/U  can  change  the  wavenumber-frequency  spectral 
level  by  3 dB) . Third,  the  speed  of  the  breathing  waves, 
and  thus  the  location  of  the  hose  resonance  is  often  only 
known  approximately.  Fourth,  it  is  difficult  to  ascertain 
that  all  significant  breathing  wave  scattering  and  other 
sources  of  extraneous  noise  have  been  eliminated  in  the 
vicinity  of  the  hydrophone  so  that  hose  resonance  is  truly 
the  dominant  noise  mechanism.  Fifth,  repeatability  of  self 
noise  measurements  between  differing  tow  platforms  and 
trials  is  about  2-3  dB  even  on  the  same  acoustic  modules. 
Finally,  agreement  between  model  predictions  and  measured 
data  on  the  order  of  3 to  5 dB  over  a frequency  band  of 
interest  is  usually  considered  to  be  quite  good.  Given 
these  uncertainities , plus  others  not  mentioned,  it  should 
not  be  surprising  that  it  is  possible  to  find  self  noise 
measurements  which  "match"  calculations  made  using  either 
mode 1 . 


The  above  discussion  indicates  that  a single 
measurement  of  the  wavenumber  spectrum  from  near  zero  to 
somewhat  beyond  the  convective  peak  should  provide 
valuable  information  on  the  wall  pressure  fluctuations 
beneath  a turbulent  boundary  layer.  Thus  a goal  of  the 
following  array  design  will  be  to  provide  such  a measurement 
over  a reasonable  frequency  range. 
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In  this  report  a very  conservative  design  philosophy 
will  be  followed.  First,  it  will  be  assumed  that  simple, 
conventional  signal  processing  techniques  will  be  used  to 
estimate  the  wavenumber-frequency  spectrum.  Second,  the 
array  design  will  be  based  on  the  model  which  places  the 
most  severe  requirements  on  the  array  performance.  In 
Figure  1.2,  the  Chase  model  has  a spectral  dynamic  range 
of  57  dB  between  the  convective  peak  and  the  low-wavenumber 
spectral  level.  This  range  will  prove  difficult  to  measure; 
thus  the  Chase  model  will  be  used  to  derive  the  array 
parameters. 


It  is  felt  that  the  array  design  developed  here 
will  be  adequate  to  measure  low-wavenumber  spectral  levels. 
However,  due  to  the  stringency  of  the  model  and  the 
assumption  of  simple  estimation  techniques  for  measuring 
the  spectral  levels,  the  resulting  array  design  will  be 
nontrivial  to  implement.  It  will  be  shown  that  approximately 
40  array  elements  are  required.  Each  element  should  be  a 
3 inch  diameter  ring  which  is  1.2  in  wide.  The  total  array 
length  should  be  about  48  inches. 

The  most  severe  requirement  is  on  the  channel-to- 
channel  tolerances.  The  rms  phase  difference  should  be  less 
than  0.29  deg.,  the  rms  gain  variations  should  be  less  than 
.06  dB,  and  rms  error  in  the  spacing  between  elements  should 
be  less  than  1.7  mils.  (Although  these  requirements  are 
decidely  nontrivial  it  should  be  remembered  that  they  represent 
rms  errors  after  calibration.  In  the  case  of  the  displacement 
errors  this  number  represents  an  rms  variation  in  the  element 
spacing.  A fixed  spacing  error  will  have  a negligible 
impact  on  the  spectrum  measurements.) 


mm, 
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Note  that  if  the  errors  are  assumed  to  be  uniformily 
distributed  between  maximum  error  limits,  the  above  rms  errors 
correspond  to  the  following  maximum  errors  — Phase:  + 0.5  deg., 
Gain:  +0.1  dB , Position:  +3.0  mils.  If  the  tolerances  are 
relaxed  by  a factor  VT,  a useful  measurement  capability  is 
still  achieved.  (See  Table  2.2) 

A future  report  will  consider  the  application  of 
high  resolution  spectral  estimation  techniques  to  reduce  th* 
required  number  of  elements. 
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2.0  ARRAY  DESIGN 


In  this  section,  the  considerations  leading  to 
the  array  configuration  described  in  the  introduction  are 
discussed.  Throughout,  it  will  be  assumed  that  the  case 
of  interest  corresponds  to  flow  along  the  axis  of  a long 
horizontal  cylinder,  and  further,  that  only  the  axisymmetric 
component  of  the  pressure  fluctuations  is  to  be  measured. 
This  allows  the  design  to  consist  of  a horizontal  array  of 
identical  flush-mounted  "ring"  transducers  as  sketched  in 
Figure  1.1.  An  estimate  of  the  axisymmetric  wavenumber- 
frequency  spectrum  can  be  obtained  by  filtering  the  output 
of  each  transducer  with  a narrowband  filter,  "beamforming" 
at  a particular  wavenumber,  and  averaging  the  resulting 
output  power.  These  operations  are  sketched  in  Figure  2.1. 


The  beamforming  operation  consists  of  multipling 
the  filter  outputs  by  a weighted  complex  phasor  and 
summing.  The  phasor  depends  on  which  wavenumber  component 
is  desired.  Thus: 

N 

e 

„ -ik  (i-l)d  \ ^ | | . .. 

= a^e  s'  ; / a | ac | = 1 (2.1) 

£=1 


where 


and 


is  the  complex  phasor  on  the  Jth  element, 

k is  the  desired  wavenumber 
s 

d is  the  element  spacing 

a?  is  a weighting  coefficient  chosen  to 

to  control  the  sidelobe  response  of  the 

* 

array 

Ng  is  the  number  of  elements  in  the  array. 


* In  general  the  weighting  coefficients  may  be  complex? 
however  for  the  current  design,  only  real  a^  are  required. 
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The  resulting  measured  wavenumber-f requency  spectrum  is 
given  by: 


PM(ks'wo)  2irlr 


J Po(k/wo) | H (k) |2|A(k-kg) |2dk 


(2.2) 


where 


Po  (k  ,u)q)  is  the  wavenumber-frequency 

spectrum  of  the  circumferentially 
averaged  pressure, 

H (k)  is  the  element  factor  representing 
the  response  of  a single  element 
to  wavenumber , k , 

A (k)  is  the  "array  factor"  given  by. 


A(k)  = 


ik(*-l)d 


(2.3) 


In  Equation  (2.2),  the  quantity  &w  is  the  white 
noise  loss  factor  given  by. 


L,  = N 
W e 


(2.4) 


The  factor  dN  represents  the  equivalent  k-space  bandwidth 
of  the  array  fictor.  That  is: 


A(k) I dk  = 


(2.5) 


-tr/d 


Note  that  the  normalization  assumed  in  (2.1)  implies  A(0)  = 1. 
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An  estimate  of  the  axisymmetric  wavenumber-frequency 
spectrum,  PQ(k,u>)  can  be  obtained  in  those  regions  of  k-space 
which  are  not  significantly  distorted  by  the  element  and  array 
factors.  Qualitatively,  the  effect  of  each  factor  in  Equation 
(2.2)  can  be  sketched  as  shown  in  Figure  2.2.  There  it  is 
shown  that  the  element  factor  is  generally  a relatively  broad 
response  in  k-space,  but  except  for  thin  elements,  significantly 
attenuates  the'  convective  component. 

The  array  factor,  on  the  other  hand,  has  a very 
narrow  response  in  k-space,  centered  on  the  desired  wavenumber 
kg.  The  approximate  width  of  the  main  response  lobe  is 
equal  to  2u  times  the  reciprocial  of  the  array  length.  By 
making  the  array  long  enough,  it  is  possible  to  obtain  high 
resolution  spectral  estimates.  If  the  wavenumber  spectrum 
were  perfectly  flat,  resolution  would  not  be  a significant 
consideration,  however,  according  to  the  model  proposed  by 
Chase  [3] , the  spectrum  changes  rapidly  with  k at  low 
wavenumbers.  If  the  spectral  level  changes  significantly 
across  the  main  response  lobe,  the  resulting  spectral  estimate 
will  be  biased  away  from  the  true  value. 

Another  feature  of  the  array  factor  which  must 
be  considered,  is  that  it  is  periodic  with  period  equal  to 
2tt  times  the  reciprocal  of  the  element  spacing.  Thus,  in 
addition  to  accepting  energy  from  the  desired  wavenumber, 
an  array  of  point  elements  accepts  energy  equally  at  the 
wavenumbers,  kg  + 2Ttn/d,  n = 1,2,...  For  a general  array, 
if  there  is  significant  energy  in  the  spectrum  at  the 
ambiguous  wavenumbers,  and  if  this  energy  has  not  been 
attenuated  by  the  element  factor,  the  resulting  measurement 
will  be  biased  by  the  ambiguities  in  the  array  factor.  In 
the  literature  on  phased  arrays,  the  array  factor  periodicities 


1 
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are  called  "grating  lobes”.  In  the  signal  processing 
literature,  the  errors  produced  by  the  periodicities  are 
called  "aliasing  errors".  Both  these  terms  will  be  used 
as  appropriate. 

Finally,  note  on  Figure  2.2,  that  between  the 
periodic  ambiguities  of  the  array  factor,  there  are  sidelobe 
responses  which  are  usually  many  dB  below  the  main  response. 
However,  if  these  sidelobes  lie  within  the  convective 
region,  they  still  may  pick  up  significant  energy  and 
consequently  bias  the  low  wavenumber  measurement. 

The  primary  problem  of  array  design  is  to  choose 
the  array  parameters  to  reduce  the  bias  errors  described 
above  to  an  acceptable  level.  The  approach  taken  in  this 
report  is  to  note  that  from  Equation  (2.2),  the  element 
factor  produces  a modified  spectrum,  P • | H. | , which  is 

measured  by  steering  the  array  factor  to  various  wavenumbers. 
The  element  pattern  produces  two  effects.  On  the  negative 
side,  it  limits  the  maximum  wavenumber  for  which  unbiased 
estimates  of  P can  be  obtained.  However,  on  the  positive 
side,  it  reduces  the  dynamic  range  of  the  spectrum  which 
is  to  be  measured  and  thus  reduces  the  requirements  on  the 
array  sidelobe  levels. 

Given  a particular  element  size,  the  product 
P • | H | is  determined,  and  then  the  problem  is  to  choose 
a high  enough  spatial  sample  rate  (small  enough  element 
spacing)  to  prevent  aliasing  errors  from  distorting  the 
spectrum.  Once  the  element  spacing  is  determined,  one 
must  choose  a long  enough  array  so  that  the  resolution 
will  be  adequate  to  prevent  biases  over  the  rapidly  varying 
portions  of  the  spectrum.  Finally,  the  array  weighting 
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function  must  be  determined  so  that  the  array  sidelobes  will 
be  low  enough  to  prevent  biases  due  to  sidelobe  leakage. 

In  general,  an  exhaustive  search  for  the  optimum 
array  characteristics  would  repeat  the  above  procedure  for 
all  possible  element  widths  until  the  best  array  was  found. 

In  this  report  the  search  will  be  conducted  in  three  stages. 

* 

First  "wafer"  elements  are  assumed.  For  such  elements  it 
is  found  that  the  maximum  frequency  for  which  useful 
measurements  can  be  made  depends  strongly  on  the  minimum 
element  separation  which  can  be  implemented.  For  a 

one-eighth  inch  separation,  the  maximum  frequency  is  160  Hz. 

** 

For  one- quarter  inch  separation,  the  maximum  is  100  Hz. 

In  the  second  stage,  a center  frequency  of  100  Hz  is 
assumed  and  the  effect  of  wider  elements  with  correspondingly 
greater  separation  is  investigated.  At  the  end  of  this 
stage,  two  nominal  array  designs  will  be  chosen.  In  the 
third  stage,  the  performance  of  the  two  nominal  designs 
will  be  investigated  for  center  frequencies  of  50  and  100  Hz. 
Also  in  this  stage,  the  effect  of  relaxed  channel-to-channel 
tolerances  is  considered.  Finally,  performance  at  other 
tow  speeds  and  frequencies  will  be  discussed. 


2 . 1 Design  for  Wafer  Elements 


Figure  2.3  plots  the  Chase  model  of  PQ  for  a three 
inch  diameter  cylinder  towed  at  15  kts.  The  results  are 


Wafer  are  so  thin  that  it  may  be  assumed  | H (k) | = 1 over 
the  region  of  interest. 

**These  limits,  of  course,  depend  on  tow  speed  and  array 
diameter,  u = 15  kts  and  2a  = 3 in  was  assumed. 
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Figure  2.3.  Wavenumber-Frequency  Spectrum  of  the  Axisymmetric 
Component  of  Pressure  for  a 3 inch  Diameter  Cylinder. 
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shown  as  a function  of  k with  frequency  of  the  narrowband 
filter,  (joq,  as  a parameter. 

Note  that  as  the  center  frequency  increases,  the 

wavenumber  corresponding  to  the  convective  peak  increases 

proportionally.  At  the  same  time  the  low  wavenumber  spectral 

content  decreases  more  rapidly  than  does  the  convective 

peak,  thus  producing  an  increased  spectral  dynamic  range. 

Both  of  these  characteristics  imply  that  the  array  design 

becomes  more  difficult  as  id  increases.  First,  as  the 

o 

convective  peak  moves  to  higher  wavenumber,  the  elements 
must  be  placed  closer  together  to  reduce  biases  produced 
by  aliasing.  Second,  the  greater  dynamic  range  requires 
better  sidelobe  control  of  the  array  response  function 
which,  in  turn,  implies  tighter  amplitude,  phase,  and 
position  tolerances  between  channels. 

Figures  2.4  and  2.5  illustrate  the  effect  of 
aliasing  for  two  cases.  These  figures  compare  the  Chase 
model  with  a numerical  evaluation  of  Equation  (2.2).  Note 
that  since  the  array  response  is  periodic.  Equation  (2.2) 
may  be  written: 


dN  yd 

5M(ks'"o>  " 2^  J PA<k' 


u>  ) | A (k-k  ) I dk 

U 5 


-TT/d 


where  PA  is  the  aliased  spectrum  given  by: 


(2.4) 


PA(k,»0) 


E 

n=-°° 


/i,  j.  2vn  , | ,,  , 2irn,  i2 

+ lH(-k  + “a-1  I 


(2.5) 


Since  only  the  effects  of  aliasing  are  of  interest,  for  this 
calculation,  a long  array  is  assumed  so  that  one  has: 
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y^lAOc-kg)  r * S(k-ks)  (2.6) 

where  6 is  the  Dirac  delta  function.  Thus  for  Figures  2.4 
and  2.5,  PM(k  ,uj  ) = P (k  ,u)  ).  Further  the  assumption  of 

MSOASO  2 * 

wafer  elements  allows  the  approximation  — | H ( k ) | = 1. 

Figure  2.4  shows  that  for  one  inch  spacing,  the 
aliased  spectrum  is  representative  of  PQ  only  in  the  region 
near  the  convective  peak.  In  this  case  the  aliased  spectrum 
PA  is  periodic  in  k with  a period  2n*(12)ft  If  the 
element  spacing  is  reduced  to  one-quarter  inch  the  period 
increases  to  2ir*(48)ft  ^ and  the  aliased  spectrum  is  a good 
representation  of  the  Chase  model.  Physically,  what  is 
happening  is  that  in  Figure  2.4,  when  the  array  is  steered 
to  low  wavenumbers,  the  first  grating  lobe  (ambiguous 
response)  of  the  array  is  folding  in  high  spectral  levels 
from  near  the  convective  peak  and  thus  masking  the  low 
wavenumber  energy.  In  Figure  2.5,  when  the  array  is  steered 
to  low  wavenumber,  the  first  ambiguous  response  is  well 
beyond  the  convective  ridge  where  the  spectral  level  is 
again  low,  and  the  folded  energy  causes  a relatively  small 
bias. 

In  evaluating  the  aliased  spectrum  (Eq.(2.5)), 
it  is  found  that  the  summation  depends  almost  entirely  on 
the  ambiguous  response  near  the  convective  peak.  Thus, 
Figure  2.3  can  be  used  to  estimate  the  required  element 
spacing  for  a three  inch  diameter  cylinder.  In  particular, 

_ 

In  the  following  figures  w denotes  the  element  width  and 
k denotes  the  effective  3 dB  array  "bandwidth"  in  k-space. 
Tnus,  D=0  is  shorthand  for  the  assumption  |H(k)|*  = 1,  and 
k = 0 is  shorthand  for  the  assumption  (2.6). 

D 
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Fi'-u’-e  2.4  . Aliasing  with  1 inch  Spacing  and  Wafer 
Elements . 
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at  low  wavenumber,  the  major  contribution  to  the  measured 
energy  will  come  from  the  ambiguous  lobe  located  at  ks+2ir/d. 

If  this  lobe  lies  too  close  to  the  convective  peak,  the  low 
wavenumber  measurement  will  be  biased  by  the  ambiguous 
response.  Quantitatively,  if  a 3 dD  bias  in  the  low  wave- 
number  measurement  is  acceptable,  the  spectral  level  at  the 
ambiguous  response  must  be  less  than  or  equal  to  the  spectral 
level  at  low  wavenumbers.  Since  the  spectrum  is  increasing 
with  wavenumber  until  the  convective  peak  is  reached,  this 
requirement  immediately  dictates  that  the  first  ambiguous 
response  must  be  well  beyond  the  convective  peak.  In 
particular,  the  element  spacing  is  chosen  so  that  the  first 
ambiguous  response  is  beyond  the  convective  peak  and  such 
that  the  spectral  level  there  is  equal  to  the  spectral  level 
at  low  wavenumber.  For  example,  at  a frequency  of  100  Hz, 
the  spectral  level  at  low  wavenumber  is  -57  dB  relative  to 
the  peak.  This  spectral  level  is  again  reached  for 
( k./2 tt ) = 47  ft  * thus  the  required  element  spacing  is  0.25  in. 

The  3 dB  criterion  represents  relatively  poor 
accuracy  for  the  spectral  estimate;  however  since  the 
discussion  in  the  introduction  indicates  that  the  current 
uncertainty  in  the  low  wavenumber  content  may  be  as  high 
as  20  dB,  it  is  felt  that  the  3 dB  criterion  is  acceptable. 

If  this  is  not  true,  the  required  spatial  sampling  rates 
will  be  significantly  increased. 

The  computation  of  the  required  sampling  rate 
based  on  the  aliasing  at  low  wavenumber  is  not  a general 
criterion  applicable  to  all  spectral  shapes.  In  general, 
the  aliasing  error  should  be  estimated  at  all  wavenumbers 
over  the  region  of  interest  and  the  spatial  sampling  rate 
should  be  chosen  to  control  the  worst  aliasing  error. 


A 
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Computations  with  the  Chase  model,  however,  indicate  that 
the  worst  aliasing  does  occur  at  low  wavenumber.  The 
above  comment  would  not  be  true  if  we  were  concerned  with 
the  measurement  of  wavenumbers  beyond  the  convective  peak. 

In  particular,  the  spectral  level  cannot  be  measured  at 
wavenumbers  greater  than  one-half  the  spatial  sampling  rate 
(Nyquist  rate)  away  from  the  convective  peak  since  the 
aliasing  error  becomes  very  large.  (See  Figure  2.4). 

Using  these  considerations,  the  required  sample 
rate  is  plotted  versus  frequency  for  a three  inch  cylinder 
in  Figure  2.6.  This  figure  demonstrates  that  there  is  a 
strong  limit  to  the  maximum  frequency  at  which  useful 
measurements  can  be  obtained.  For  example,  if  it  is  assumed 
that  an  element  spacing  of  one-eighth  inch  represents  the 
minimum  which  could  be  achieved,  then  the  maximum  frequency 
is  limited  to  about  160  Hz.  For  one-quarter  inch  spacing 
the  maximum  frequency  is  about  100  Hz. 

The  next  step  in  determining  the  array  parameters 
is  to  compute  the  required  resolution  in  k-space  which  will 
be  related,  later,  to  the  required  array  length.  To  make 
this  computation,  the  main  response  lobe  of  the  array  factor 
is  approximated  by  the  Gaussian  function  over  a limited 
region  of  k-space  near  the  desired  wavenumber  k^.  That  is, 

' lk-ksl  i 3°k 

(2.7) 

where  o.  = (8  ln2)  * ^k_ 

K D 

and  kn  is  the  3dB  "bandwidth"  in  k-space. 

D 
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Using  this  approximation,  the  measured  spectral  level  can 
be  computed  using  Equation  (2.2)  and  a numerical  quadrature 
formula  to  evaluate  the  integral.  Since  the  quantity  being 
investigated  here  is  the  bias  produced  by  the  finite 
resolution  of  the  array  factor,  the  integration  is  restricted 
to  the  main  response  lobe  near  the  desired  wavenumber  ks- 
The  effect  of  sidelobe  response  will  be  considered  separately 
below.  The  Gaussian  approximation  to  the  main  response 
lobe  is  quite  good  independent  of  which  array  weighting 
coefficients  are  used.  The  major  effect  of  the  weighting 
coefficients  is  to  control  the  sidelobe  level  at  the  expense 
of  increased  k-space  bandwidth.  Thus  different  array 
coefficients  require  different  array  lengths  for  the  same 
k . As  long  as  the  results  are  given  in  terms  of  resolution, 
they  are  quite  general. 

The  result  of  a numerical  evaluation  of  Equation  (2.2) 

* 

using  approximation  (2.7)  for  A(k)  is  shown  in  Figure  2.7. 

For  that  computation  the  elements  were  assumed  to  be  0.25  in. 
apart,  and  the  array  bandwidth  was  taken  as  ( k^/ 2tt ) = 1.0  ft 
Again,  the  most  significant  bias  occurs  near  zero  wavenumber. 

If  the  bias  is  plotted  versus  array  bandwidth,  for  kg  = 0, 

Figure  2.8  results.  From  that  figure,  it  can  be  seen  that 
a resolution  of  (kn/2ii)  = 0.65  ft  1 is  required  to  maintain 
less  than  3 dB  bias  relative  to  the  Chase  model.  The 
required  array  length  will  depend  on  the  particular  weighting 
coefficients,  and  these,  in  turn,  depend  on  the  required 
sidelobe  level. 

Figure  2.9  can  be  used  to  estimate  the  required 
sidelobe  level  of  the  array  factor.  As  can  be  seen. from 


For  this  and  subsequent  computations,  the  "wafer"  elements 
will  be  assumed  to  be  one  quarter  inch  wide  (w-.25in)  as  this 
produces  negligible  distortion  of  the  spectral  measurement 
over  the  regions  of  interest. 
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Figure  2.7.  Effect  of  k-space  Bandwidth  on  the  Measured 
Spectrum,  P_. . 
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that  figure,  this  level  is  a function  of  the  dynamic  range 
of  the  spectrum,  DR,  its  k-space  bandwidth  near  the  convective 
peak,  B,  and  the  array  bandwidth  kg.  Specifically,  the  ratio 
of  the  leakage  through  the  array  sidelobes,  P , to  the 

b L 

spectral  level  measured  at  low  wavenumber  is  given  approximately 
by:  * 


10 


SL  - DP.  + 10  log  (k_./B) 


(2.8) 


where  SL  is  the  rms  sidelobe  level  near  the  convective  peak. 

If  the  sidelobe  leakage  is  to  be  less  than  the  measured 
value,  the  sidelobe  level  must  satisfy: 

SL  > DR  + 10  log (H/kg)  (2.9) 

For  a 3 inch  cylinder  towed  at  15  kts  and  a center 
frequency  of  100  Hz,  the  dynamic  range  is  57  dB  and  the 
bandwidth  of  the  convective  peaks  is  ( B/2 tt ) = 1.2  ft 
Substituting  these  values  into  (2.9)  yields  a sidelobe 
requirement  of  about  60  dB. 

Figure  2.10  illustrates  the  effect  of  various 
sidelobe  levels  using  a numerical  evaluation  of  Equation  (2.2). 
For  the  numerical  evaluation,  a constant  sidelobe  level 
was  assumed.  Figure  2.10  is  consistent  with  the  approximation 
(2.9)  as  it  indicates  that  a sidelobe  level  of  60  dB  or 

it 

To  be  specific,  the  dynamic  range  is  defined  as  DR  = 

10  log  [Po(uj/u  ,u>) /P  (0  ,io)  ] and  B is  the  bandwidth  between 
the  3 dB  points  around  the  convective  peak.  The  rms 
sidelobe  level  is  defined  as  SL  * - 10  log  [ (d/2n)  / | A (k)  | Zdk] , 
where  the  region  of  integration  is  confined  to  the  sidelobe 
region  of  A (k) . 
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Figure  2.10.  Effect  of  Sidelobes  on  the  Measured  Spectrum. 
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greater  is  required  to  prevent  sidelobe  leakage  from  the 
convective  peak  from  masking  the  low  wavenumber  content. 

This  represents  a severe  sidelobe  requirement. 
Although  in  theory,  Dolph-Chebyshev  weighting  functions 
[7],  [8],  achieve  such  low  levels,  tight  control  on  channel- 

to  channel  amplitude,  phase,  and  position  errors  is  required. 
Table  2.1  summarizes  the  important  characteristics  of  the 
Dolph-Chebyshev  weights.  This  class  of  weighting  functions 
should  be  useful  for  the  current  application  as  they  achieve 
a given  peak  sidelobe  level  with  minimal  broadening  in  the 
main  response  lobe.  Table  2.1  lists  the  "rms"  sidelobe 
level,  the  3 dB  bandwidth  and  the  white  noise  loss  factor 
for  Dolph-Chebyshev  weights  as  a function  of  the  peak 
sidelobe  level.  From  that  table  it  appears  that  the  array 
length  required  to  achieve  60  dB  rms  sidelobes  and  a k- 
space  bandwidth  of  2v • (.65)  ft  ^ is  2.2  ft.  With  an 
element  spacing  of  0.25  inches,  this  corresponds  to  105 
elements . 


As  noted  earlier  the  above  parameters  represent 
a severe  requirement  on  the  amplitude,  phase,  and  position 
tolerances  between  the  channels.  It  is  shown  in  Appendix  B 
that  the  effect  of  independent  amplitude  and  phase  and 
position  errors  on  the  array  response  is  given  by: 


E|B(k,kg) | 


|Ao(k-ks)|2(1-o2-k2o2)  . 


N 


(2.10) 


where 


2 

E | B (k ,k  ) | is  the  statistical  expectation  of 
the  array  response  (squared)  which 
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TABLE  2.1 

ARRAY  FACTOR  CHARACTERISTICS 
FOR  DOLPH-CHF.BYSHF.V  WEIGHTING 


COEFFICIENTS 


Peak.  Sidelobe 
Level 
(dB) 


rms 

Sidelobe  Level 
(dB) 

(approx) 


Normalized 
White  Noise 
Loss 
Factor 


25 


28 


1.09 


30 


33 


1.14 


35 

40 

50 

60 

70 

80 


38 
43 
53 
6 3 
73 
83 


1.21 

1.29 

1.39 

1.51 

1.62 

1.73 


Normalized 

3-dB 

Bandwidth 


1.01 

1.07 

1.13 

1.20 

1.33 

1.44 

1.55 

1.65 
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has  been  perturbed  by  the  amplitude 

phase  and  position  errors, 

is  the  unperturbed  array  factor, 

(It  is  assumed  that  A^(k)  has  been 

normalized  so  that  A (0)  = 1) 

o 

is  the  white  noise  loss  factor  for 
the  particular  weights  used  (column 
3 in  Table  2.1) 
is  the  number  of  elements, 
is  the  rms  phase  error, 
is  the  rms  amplitude  error  as  a 
ratio  (i.e.,  t^ffia^J/a.)  , 
is  the  rms  positional  error. 

The  60  dB  sidelobe  requirement  must  be  satisfied 
at  the  convective  peak  with  the  array  steered  to  zero,  thus 
one  must  have  F.  I Blw/15^,0)  ( “ < 10  . If  the  errors  are 

distributed  "equally"  between  the  three  sources,  one  must 
have : 

— o_  =0=0=  0.0049  (2.11) 

u ts  c 

c 

That  is,  the  phase-error  must  be  less  than  .28  deg.  rms, 

* 

the  amplitude  error  must  be  less  than  0.042  dB  rms  and  the 
position  error  must  be  less  than  1.6  x 10  ? in. 


A (k-k  ) 
o s' 


W 


N 


and 


% 


It  is  generally  useful  to  give  amplitude  tolerances  in 
dB  as  these  tolerances  are  directly  relatable  to  tolerances 
on  the  preamp  and  amplifier  gains  on  each  hydrophone  channel. 
To  convert  to  variance  of  the  dB  value  (say  of  the  amplifier 
gain)  note  that  for  small  errors,  gain  variations  are  related 
to  fractional  amplitude  variations  as: 


AC 


a ( a { 

20  log  — = 8.68  In  — 


8.68  ln(l+t~f)  ~ 8.68t, 
8 . 6 8 n 

c 


thus 
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In  summary,  the  design  based  on  wafer  elements 

requires  105  one-quarter  inch  elements  at  the  same  spacing. 

The  array  length  is  about  2.2  ft.  Sidelobe  levels  of  60  dB 

are  required  which  imply  channel-to-channel  tolerances  of  — 

o_  = .042  dB , o.  = .29  deq . , o,  = 1.6  mils.  These  tolerances 
'J  T P 

will  be  difficult  to  achieve.  It  should  be  remembered 
however,  that  they  represent  residual  variations  after 
calibration.  A fixed  amplitude,  phase  or  spacing  error 
will  have  a negligible  effect  on  the  spectral  measurement. 

If  the  errors  are  uniformly  distributed  between  maximum 
limits  the  above  rms  errors  correspond  to  | AG  | <_  .07  dB , 

| A<f  | <^0.5  deq.,  | A 3 1 ^ 2.8  mils.  The  required  number  of 
channels  is  larqe,  and  the  k-space  bandwidth  of  (kQ/2Ti)  = 0.65  ft 
is  marginal  as  will  be  seen  later. 

2 . 2 Effect  of  Element  Length 

The  results  of  Section  2.1  indicate  that  an  array 
design  with  wafer  elements  will  require  both  tight  channel- 
to-channel  tolerances  and  a large  number  of  elements  to 
achieve  a qood  low-wavenumber  measurement  capability.  In 
this  section,  designs  based  on  wider  elements  will  be 
considered.  Wider  elements  have  a narrower  response  in 
k-space,  which  has  the  effect  of  reducing  the  required 
spectral  dynamic  range,  at  the  expense  of  reducing  the 
maximum  k to  which  unbiased  measurements  can  be  made. 

In  [1] , two  approaches  to  rejecting  the  convective 
peak  were  considered.  The  first  approach  was  based  on 
placing  the  peak  far  out  on  the  sidelobes  of  the  element 
pattern.  This  approach  was  rejected  as  it  leads  to  very 
long  elements  (5-20  ft). 


2-26 


The  second  approach  was  based  on  choosing  the 
element  width  to  place  the  first  null  of  the  element 
pattern  near  the  convective  peak.  The  objective  was  to 
reduce  the  spectral  dynamic  range  and  thus  ease  the  channel- 
to-channel  tolerance  requirements.  The  argument,  unfortunately, 
is  fallacious. 

As  an  example,  it  was  found  that  at  U = 15  kt, 

2a  = 3 in,  and  (lo/2tt)  = 100  sec  a two  inch  element  was 
required  to  place  the  first  null  near  the  convective  peak. 
However,  to  prevent  significant  aliasing  the  modified 
spectrum  had  to  be  over  sampled  at  one  to  two  samples  per 
inch.  It  was  proposed  that  the  over  sampling  could  be 
accomplished  by  summing  overlapping  groups  of  smaller 
elements  thus  giving  an  effective  element  size  of  two  inches 
with  samples  spaced  less  than  two  inches  apart.  For 
example,  if  one-half  inch  physical  elements  are  used, 
overlapping  groups  of  four  can  provide  an  effective  element 
length  of  two  inches  with  samples  every  one-half  inch. 

The  fallacy  is  that  now,  the  array  consists  of 
one-half  inch  elements  not  two  inch  elements.  The  spectral 
dynamic  range  is  that  of  half-inch  elements.  The  correct 
way  of  viewing  the  problem  is  to  note  that,  the  process  of 
summing  overlapping  groups  of  elements  and  then  applying 
weighting  factors  to  the  groups  is  equivalent  to  applying 
a modified  set  of  weights  to  the  original  element  outputs. 

The  only  virtue  of  the  modified  weights  is  that  they  have 
a very  low  sidelobe  level  in  the  neighborhood  of  the 
convective  peak.  However,  their  tolerance  to  amplitude, 
phase,  and  position  errors  is  no  better  than  say  Dolph- 
Chebyshev  weights.  In  particular,  note  that  from  Equation  2.10, 
the  major  dependence  of  the  sidelobe  error  on  the  weights  is 
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1 


in  the  term  ■ . This  quantity  does  not  change  greatly 

between  different  sets  of  weights;  thus  the  bias  in  the 
array  factor  will  be  about  the  same  for  the  modified  weights 
as  for  an  equivalent  set  of  Dolph-Chebyshev  weights.  The 
tolerance  requirements  nre  also  the  same.  In  the  following 
discussion,  it  w: l 1 be  assumed  that  the  element  width  and 
spacing  are  equal. 

The  first  step  in  estimating  the  effect  of  the 
element  width  is  to  compute  the  difference  between  the  Chase 
model  and  the  aliased  spectrum  (as  modified  by  the  element 
factor,  Eq.  2.5).  Figure  2.11  illustrates  the  results  of 
this  computation.  There  the  region  in  k-space  over  which 
the  aliased  spectrum  is  not  more  than  3 dB  greater  than 
the  Chase  model  is  plotted  versus  element  width.  For  small 
elements,  the  region  of  unbiased  measurement  is  | lc/2  tt  | <,  1/2W. 
Note  that  as  the  element  width  and  spacing  increases,  the 
region  of  valid  measurements  is  correspondingly  reduced. 

Since  the  most  of  the  spectral  energy  is  concentrated  near 
the  convective  peak,  the  measurement  interval  tends  to  be 
centered  around  this  peak.  For  one-inch  elements,  the 
convective  peal;  is  reduced  about  3 dB.  For  larger  elements, 
the  aliasing  error  increases  rapidly  and  the  measurement 
interval  is  rapidly  limited  to  the  acoustic  region.  The 
maximum  element  size  which  yields  a measurement  interval 
beyond  the  acoustic  region  is  about  1.4  in. 

Figure  2.12  shows  the  spectral  dynamic  range  versus 
element  size.  Note  that  this  range  does  not  decrease 
significantly  until  the  element  size  exceeds  1.4  in.  That 
is,  if  the  element  length  is  chosen  large  enough  to 
significantly  reduce  the  spectral  peak,  the  corresponding 
aliasing  errors  have  reduced  the  measurement  interval  to 
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Figure  2.12.  Spectral  Dynamic  Range  Versus  Element  Size. 
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the  acoustic  region. 

Figure  2.13  shows  the  number  of  channels  required 
to  achieve  (k_/2TT)  = 0.65  ft  ^ and  Figure  2.14  shows  the 
required  channel--to-channel  tolerances,  Dolph-Chebyshev 
weights  are  assumed.  The  least  severe  tolerances  for  the 
least  number  of  elements  occurs  at  w = 1.4  in.  This 
unfortunately  is  a spurious  optimum  since  its  usefulness 
depends  on  carefully  tailoring  the  element  pattern  to  the 
spectrum  to  be  measured.  For  the  remaining  calculations 
a somewhat  more  conservative  maximum  element  width  of  1.2  in. 
is  assumed. 


At  the  other  end,  the  channel  tolerances  can  be 
eased  by  using  more  elements.  However , the  tolerances 
increase  as  1/  ; thus  a large  number  of  channels  is 

required  to  significantly  relax  the  tolerances.  Note  that 
decreasing  the  element  size  from  one  inch  to  one-half  inch 
has  a negligible  effect  on  the  required  tolerances.  This 
is  due  to  a 2 dB  increase  in  the  spectral  dynamic  range 
which  nearly  compensates  the  increased  number  of  elements. 

Figures  2.13  and  2.14  do  not  indicate  the  true 

tradeoffs  between  element  width,  number  of  elements  and 

channel-to-channel  tolerances.  In  particular,  those  figures 

were  constructed  for  a constant  k-space  bandwidth  of 

(k„/2iT)  = 0.65  ft  "*■.  Recall  that  this  is  the  maximum 
B 

allowed  k-space  bandwidth.  There  is  no  penalty  incurred 
for  decreasing  this  bandwidth  by  using  a longer  array.  In 
fact,  from  Figure  2.8,  decreasing  the  bandwidth  will  reduce 
the  bias  due  to  nonzero  mainlobe  width  and  improve  array 
performance.  Thus  if  the  channel-to-channel  tolerances 
must  be  relaxed  by  using  more  channels,  it  is  more  advantageous 


M 
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Figure  2.13.  Number  of  Channels  Required  to  Achieve 
(k_/2ir)  = 0.65  ft"1. 
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Figure  2.14.  Required  Channel-to-Channe)  Tolerances  versus 
Element  Size. 


to  use  the  largest  possible  element  size  and  increase  the 
array  length  accordingly  rather  than  use  smaller  elements. 

This  is  demonstrated  in  Figure  2.15,  where  the  phase 
tolerance  is  plotted  versus  the  number  of  elements  with 
element  width  as  a parameter. 

Figure  2.15  gives  the  relation  between  the  required 
number  of  channels  and  the  required  tolerances.  If  either 
limit  is  known  the  other  may  be  estimated  from  that  figure. 
Currently  neither  limit  is  known.  Flush  mounted  elements 
of  the  type  indicated  in  Figure  1.1  have  never  been  built 
and  calibrated  in  situ,  therefore  there  is  no  directly 
applicable  previous  experience  which  can  be  used  to  estimate 
achievable  tolerance  limits.  Secondly,  at  this  preliminary 
design  stage,  it  is  not  known  whether  existing  tow  cable 
and  telemetry  systems  will  be  used  or  new  ones  will  be 
built.  Thus  the  number  of  channels  available  is  not  known. 

The  approach  taken  here  will  be  to  "ballpark" 
the  achievable  tolerances  with  optimistic  estimates.  In 
particular,  it  will  be  assumed  that  the  maximum  gain 
variation  can  be  held  to  + 0.1  dB,  the  maximum  phase  variation 
is  + 0.5  deg.  and  the  maximum  relative  position  error  is 
+ 3 mils.  If  the  errors  are  uniformly  distributed 
between  their  maximums,  the  corresponding  rms  errors  are: 
aQ  = .058  dB,  = .29  deg.,  and  a ^ = 1.7  mils.  These 
tolerances  are  optimistic,  the  gain  and  phase  tolerances 
are  two  to  three  times  tighter  than  typical  specifications 
for  acoustic  modules.  The  justification  for  assuming  these 
low  values  is  that  they  are  intended  to  represent  tolerances 
after  calibration  corrections  have  been  applied. 

From  Figure  2.15,  the  above  tolerances  imply 
35  to  40  elements  will  be  required  to  make  low  wavenumber 
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measurements.  This  fact  also  influenced  the  choice  of 
the  assumed  tolerance  limits  as  it  is  unlikely  that  the 
number  of  channels  available  will  greatly  exceed  40. 

For  comparison,  the  performance  of  a design  based  on  20 
elements  with  tolerances  which  are  smaller  by  V2  will  also 
be  considered. 

Figures  2.16  through  2.23  plot  the  measured 

spectrum  as  computed  from  a numerical  evaluation  of  Equation  2.2 

in  comparison  with  the  Chase  model.  Evaluations  were  made 

for  two  arrays  at  two  different  frequencies  — 50  Hz  and 

100  Hz.  Additionally,  computations  were  made  for  sidelobe 

levels  corresponding  to  the  "desired"  tolerances  and  for 

tolerances  which  have  been  increased  (relaxed)  by  VI.  m 

those  figures,  a is  the  common  tolerance:  a = = 0£ 

= (o)/U  )oD.  Note  also  that  the  abscissa  is  a "lin-log" 

® p .. 

scale.  That  is,  the  scale  is  proportional  to  (sgn  k)log|k| 

for  | k/2 tt  | > .01  and  is  linear  for  | k/2 7T  | £ 0.01,  where 

"sgn"  is  the  signum  (or  sign)  function. 

First,  it  should  be  noted  that  the  "prefered" 
array  at  its  design  frequency  of  100  Hz  has  a bias  of  about 
6 dB  in  the  low-wavenumber  measurement.  (See  Figure  2.16) 

Since  the  array  was  designed  to  achieve  a nominal  3 dB 
accuracy  or  better  it  is  reasonable  to  ask  why  the  bias 
is  so  large.  The  answer  is  that  the  array  was  designed 
so  that  each  component  of  the  total  bias  was  3 dB  or  less. 

The  total  error  is  greater.  In  particular,  it  can  be  shown 
that  at  zero  wavenumber  the  bias  due  to  aliasing  is  0 dB , 
the  bias  due  to  the  finite  k-space  bandwidth  is  1.7  dB 
(See  Figure  2.8)  for  (kfi/2Tr)  = .36  ft  1 and  the  bias  due 
the  sidelobe  leakage  is  3.5  dB  (Eq.  2.8  and  2.10).  The 
"power  sum"  of  these  biases  is  5.7  dB  as  seen  in  Figure  2.16. 
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Figure  2.16.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Large  Array,  Tight  Tolerance,  100  Hz). 
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Figure  2.17.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Large  Array,  Tight  Tolerances,  50  Hz). 
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Figure  2.18.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Large  Array,  Relaxed  Tolerances,  100  Hz). 
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Figure  2.19.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Large  Array,  Relaxed  Tolerances,  50  Hz). 
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Figure  2.20.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Small  Array,  Very  Tight  Tolerances,  100  Hz). 
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Figure  2.21.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Small  Array,  Very  Tight  Tolerances,  50  Rzi. 
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Figure  2.22.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Small  Array,  Very  Tight  Tolerances,  100  Hz). 
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Figure  2.23.  Comparison  of  the  Measured  Spectrum  with  the 
Chase  Model  (Small  Array,  Very  Tight  Tolerances,  50  Hz). 
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At  50  Hz,  (Figure  2.17)  the  preferred  array 
performs  well.  It  provides  a measurement  within  3 dB  of 
the  Chase  model  over  the  range  -2ft  * <_  (k/2n)  4 ft 

This  range  includes  the  convective  peak  at  (k/2ir)  = 2.9ft 
Even  when  the  tolerances  are  relaxed  by  V2  , a useful 
measurement  range  is  achieved  --  0.2  ft  * ^ ( k / 2 tt ) _<  2 ft  ^ . 
The  measured  value  at  the  convective  peak  is  about  4.5  dB 
below  the  Chase  model.  The  good  performance  at  50  Hz 
can  be  attributed  to  the  fact  that  the  spectral  dynamic 
range  is  6 to  7 dB  less  at  this  frequency  than  at  100  Hz. 

As  can  be  seen  from  Figure  2.18,  the  performance 
at  100  Hz  is  measurably  degraded  by  the  relaxed  tolerances. 
The  bias  at  low  wavenumber  increased  from  5.7  dB  to  7.5  dB. 

Figures  2.20  through  2.23  indicate  the  performance 
achieved  by  the  shorter  array.  Due  to  its  greater  k-space 
bandwidth  the  biases  have  increased  even  though  the 
tolerances  have  been  decreased  to  maintain  the  same  sidelobe 
level.  The  k-space  region  for  which  the  bias  in  the 
measured  spectrum  is  3 dB  or  less  is  reduced  to  .6  ft.  ^ 

(k/2n)  < 4 ft  1 for  a center  frequency  of  100  Hz.  The 
maximum  bias  at  low  wavenumber  is  6.5  dB.  These  numbers 
indicate  performance  comparable  to  the  40  element  array, 
but  with  greatly  relaxed  tolerances  on  the  larger  array. 

(i.e.  with  the  olerances  on  the  longer  array,  greater  by 
a factor  of  two) . 

At  50  Hz  the  shorter  array  performs  significantly 
poorer  than  the  longer  array.  The  measurement  interval 
has  been  reduced  to  a relatively  narrow  region  around  the 
convective  peak.  At  50  Hz,  the  Chase  model  predicts 
a more  narrow  spectrum  and  the  wider  k-space  bandwidth 
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of  the  20  element  array  is  insufficient  to  reproduce  it 
adequately. 

The  reduced  element  tolerances  do  not  greatly 
further  degrade  the  shorter  array.  This  due,  primairly, 
to  the  fact  that  the  shorter  array  is  already  significantly 
limited  by  its  k-space  bandwidth. 

Table  2.2  summarizes  Figures  2.16  through  2.23. 

That  table  indicates  why  the  40  element  array  is  preferred. 

At  50  Hz  with  the  tighter  tolerances  it  provides  a very  good 
measurement.  It  has  the  capability  to  reproduce  the  Chase 
model  (within  3 dB)  from  negative  wavenumbers  to  wavenumbers 
beyond  the  convective  peak.  At  100  Hz,  the  measurement 
interval  is  still  reasonably  good.  The  low  wavenumber 
limit  (k/2Tr  = .4  ft  ^)  corresponds  to  wavespeeds  on  the 
order  of  250  ft/sec.  This  limit  is  comparable  to  the 
breathing  wavespeed  for  many  types  of  hose  construction; 
thus  the  measurement  provides  useful  information.  The 
upper  limit  ( (k/2Tr)  = 4 ft  1)  is  somewhat  below  the  convective 
peak.  If  the  tolerances  are  relaxed  by  & (aG  = .082  dB, 

= .409  deg,  Og  = 2.45  mils,  SL  = 51.5  dB)  a useful 
measurement  range  is  still  achieved,  particularly  at  the 
lower  frequency. 

The  20  element  array  with  very  tight  tolerances 
(aG  = .041  dB,  = .289  deg,  Ug  = 1.73  mils)  provides  a 
moderate  capability  at  100  Hz,  but  is  quite  limited  at 
50  Hz  due  tc  its  relatively  wide  k-space  bandwidth. 
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2 . 3 Capability  at  other  Tow  Speeds  and  Frequencies 

The  last  two  columns  of  Table  2.2  indicate  the 

Strouhal  number,  N , and  normalized  wavenumbers,  ka, 

s 

corresponding  to  the  frequency  and  measurement  intervals 
respectively.  Since  the  Chase  model,  as  presented  in 
Appendix  A,  depends  only  on  Strouhal  number  and  normalized 
wavenumber,  the  prefered  array  should  provide  a good 
measurement  capability  for  N <,  3,  and  0 < ka  <,  .08.  (Of 
course  if  the  array  diameter  is  changed,  the  element  size 
should  be  scaled  proportionally) . The  above  limits  are 
reasonable  engineering  "guesstimates",  however  the  scaling 
of  the  spectrum  given  in  Appendix  A should  not  be  accepted 
without  reservation.  In  particular,  a dimensional  analysis 
indicates  that  the  spectrum  should  depend  on  other  dimensionless 
parameters  such  as  (v*/U) , (6/a),  etc.  The  model  presented 

in  Appendix  A makes  the  implicit  assumption  that  for  typical 
towspeeds , array  size  and  construction,  etc.  these 
parameters  are  nominally  constant.  This  is  probably  incorrect, 
in  general,  but  is  a reasonable  engineering  approximation 
for  three  inch  cylinders  at  typical  tow  speeds. 
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APPENDIX  A:  SPECTRAL  NOTATION 

Let  p(c,0,x)  denote  the  pressure  fluctuations  on  the 
surface  of  the  cylinder  at  distance  z along  the  cylinder, 
angle  0 around  the  cylinder  and  time  t.  The  pressure  is 
assumed  to  be  stationary  in  all  three  coordinates  so  that 
the  origin  of  the  coordinate  system  is  not  important. 

Let  W(Z,6,t)  denote  the  space-time  correlation  function: 


K U , 0 , t ) = F(p(;1,01,t1)P(?1+C,01+0,'[1+t)  ) 


(A-l) 


where  "E"  denotes  statistical  expectation.  Since  W is 
periodic  in  0 with  period  2n,  it  has  a mixed  Fourier-transform 
Fourier-series  representation: 
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with  the  reciprocal  relation: 
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In  the  above  relations,  P (k,oj)  is  the  roth  angular 

m 

harmonic  spectral  density  in  streamwise  wavenumber  and 
frequency.  The  total  mean  square  pressure  fluctuation 
is  given  by: 
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The  quantity  Po(k,w)  is  the  wavenumber  frequency  spectrum 
of  the  axisymmetric  component  of  pressure. 

Here  the  narrowband  coherence  function  is  defined  as: 
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and  therefore  is  related  to  the  angular-harmonic  spectral 
densities  by: 
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In  this  notation,  the  frequency  spectral  density  of 
the  directly  transmitted  TBL  pressure  in  the  center  of  a 
fluid  filled  hose  is  given  approximately  by: 
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where  G (k)  is  the  hose  transfer  function  from  axisymmetric 
o 

pressure  to  internal  pressure.  For  a lightly  damped  hose, 
with  a breathing  speed  much  greater  than  the  convection 
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velocity.  This  is  approximated  as: 
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where  c is  the  breathincr  wavespeed  and  C is  the  hose 
o 


damping  facto*  The  one-sided  spectral  level  (per  Hz) 
is  given  as . 


LT(a>)  = 10  log  M ?T  ( ■>  1 
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Chase  [2]  gives  the  following  expression  for  Pq: 
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c t = 0.063 


v*  = . 04U 


K2  = k2  + 
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U = 0.68U 
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A = 1.08a 
h = 3.7 


In  general,  the  various  dimensionless  parameters, 
cfc,  (v*/U) , (A/a),  h,  (Uc/U) , etc.  are  thought  to  depend 
on  (6/a) , where  5 is  the  boundary  layer  thickness.  The 


A- 4 


ratio  ( /a) , in  turn,  will  depend  on  downstream  location. 
However , Chase  does  not  give  functional  relations  for  these 
parameters.  Here  the'.'  are  assumed  to  be  constant  over  the 
applicable  range  of  tow  speeds  and  streamwise  position. 

Kith  this  vor'j  crude  assumption,  one  may  write: 


2 3 ~> 


and 


Po(k'  1 ^t^v*a"(Ka)“(Kta) 


(Ka)2  = (ka)2  a , * 


(A-ll) 


(K. 


/ 

1 l'  -Vi 

V J 

\ \ 

, 2 

((  UAhv  J 

UvJ 

')  * 

(ka)  + .85 

and  the  spectrum  scales  according  to  the  Strouhal  number 
Ng  = (oa/U)  and  normalized  wavenumber,  ka. 

Gardner  [3]  used  a different  normalization  in  developing 
his  model  for  Pq.  If  this  author  has  correctly  interpreted 
the  normalization  used  by  Gardner,  his  model  may  be  written: 


P_(k,io)  = £ 

8 


(u)  f i-e~TTOt2N8 

,2  [^«2Ns 


(A-12) 


where 


2a,  N 
1 s 


f^f+(ka  - iK)1’ 


In  the  above  expression  Ng  is  the  Strouhal  number  (tua/U)  , 
and  ct^  and  oi^are  constants  given  by  --  = 0.1,  a 2 ~ 1.7. 

Again  making  the  crude  assumption  that  , ci2 , (v  /U) , 
and  Uc/U  are  constant,  the  Gardner  model  also  scales 
according  to  and  ka. 

As  a check  on  the  normalization,  it  is  noted  that 
Equation  (A-9)  agrees  with  that  given  by  Gardner  [3]  for 
the  directly  transmitted  TBL  pressure. 

For  purposes  of  scaling,  the  plots  presented  earlier, 
were  normalized  by  the  maximum  over  frequency  and  wavenumber 
of  the  Chase  model.  This  maximum  is  given  approximately 
by : 

P__„  = max  max  P (k,w) 

ma  x o 

— t»<k<ot’  — °°<  to  <«> 


2 3 2 

0. 0172p  v*a 


(A-13) 


B-l 


APPENDIX  B:  THE  EFFECT  OF  AMPLITUDE,  PHASE  AND  POSITION 
ERRORS  ON  THE  ARRAY  FACTOR 


The  array  response  to  wavenumber,  k when  steered  to 

wavenumber  k can  be  written  as: 
s 


B(k.ks)  =■  £ vik[(«-l)d+8l] 


(B-l) 


where 


g.  is  the  position  error  of  the  fcT"  element. 


Now  from  Equation  2.1: 


V*  = a*e 


-ik  ( Jt  — 1 ) d 


(B-2) 


so  that 


B(k,k3)  - £ ateil(k-kS>“-1,d  + k6li 


(B-3) 


Note  that  B cannot  be  written  as  a function  only  of  the 

difference  k-k  . 

s 

Now  if  amplitude  and  phase  errors  are  present,  a^ 
can  be  written: 


az  = a°  ( 1+e  e1^#. 


(B-4) 


where 


a°  is  the  true  (error  free)  amplitude, 
e.  is  the  fractional  amplitude  error 


E- 1 


and 

$ . is  th  phase  error. 

Substituting  (H  -i ) in'o  U 3),  and  -guaring  the  magnitude 

yields : 


| B ( k , k s ) |2 


I ( ( -iii ) i.l 


d+f*> 


(i+e  )ei{W 

m 


(B-5) 


where  v is  given  by: 


* 16 


kj? 


e 


(B-6) 


Assuming  zero  mean  errors,  independence  between  amplitude, 
phase  and  position  errors,  independence  between  channels 
and  small  errors  yields: 


E ( 1+e  p ) + m5  = E [ (1+e., ) (l+tm)  ] Eej  (','lT','m) 


m 


E|1  + * j2(W2/21 


1 + lW2  - n-W’J 


(B-7) 


Thus  the  expected  value  of  the  magnitude  of  the  array 
response  is  given  v: 
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E]B(k,ks) 


o o ik(£-m)d 
a n a_e 


£ ,m 


-c 


+ 

l 

‘ O o ik ( e-m) d 

. a n a e 

st  m 


= lA(k-ks)  |2(l-a2-k2a2)  + (a2+a2+k2a2)  £(aJ)2 


e-  ~ 


(B-8) 


Recall  that  the  white  noise  loss  factor  is  given  by: 


Ne  E (a°> 


o,  2 


V = 


(E  a°) 


o>  2 


N 


o,  2 


(B-9 ) 


since  it  has  been  assumed  that 


A(0)  = 


= 1 


(B-10) 


Thus  the  average  power  response  can  be  written: 


E|B(k,ks)|2  = |A(k-ks) |2(l-o2-k2a2) 


+ {°l+%+k2^  jT 

e 


(B-ll) 
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LIST  OF  SYMBOLS 


B(k,kg) 


oo(k, 


cylinder  radius 

weighting  coefficient  chosen  to  control 
the  side lobe  response  of  the  array 
errorless  weighting  coefficient 
"array  factor"  (Eq.  2.J) 

k-space  bandwidth  near  the  convective  peak 
array  response  at  wavenumber  k when  steered 
to  wavenumber  k 

s 

acoustic  propagation  velocity 
breathing  (bulge)  wavespeed 
dimensionless  constant  in  the  Chase  model 
of  PQ(k,toi) 

separation  between  elements 

spectral  dynamic  range  (in  dB) 

statistical  expectation  operator 

reciprocal  period  (frequency) 

hose  transfer  function  from  axi symmetric 

pressure  to  internal  pressure 

element  factor  representing  the  response  of 

a single  element  to  wavenumber,  k 

wavenumber  (2  it  times  reciprocal  wave  length) 

k-space  3 dB  bandwidth  of  the  array  factor 

wavenumber  at  convective  peak  (k  = co/u  ) 

c c 

wavenumber  to  which  the  array  is  steered  (Eq.  2.1] 

white  noise  loss  factor  (Eq.  2.4) 

one  sided  spectrum  level  of  the  directly 

transmitted  TBL  pressure  (Eq.  A-9) 

number  of  array  elements 

Strouhal  number,  (ma/u) 


S-2 


t) 

PA 

’’m 

r (k.ui) 

i’T  («•-') 


Q(u) 


SL 


U 

U 

c 


w 


w ( c , e , r ) 


A 

AG 

AH 

\>t> 

f « 


r 


o 


(0 


point  pressure 
sl.iased  spectrum  ( Eq . 2.t>) 
niLMburod  wavenumber  -frequency  spectrum 
wavenumber- frequency  spectrum  for  the 
a s i symme t i i c component  of  wall  pressure 
frequency  spectral  density  of  the  directly  TBL 
pressure  (Eq.  A- 7) 

measured  frequency  spectral  density  of  the 

point  surface  pressure  under  a TBL 

array  side  lobe  level  (in  dB) 

free  stream  speed 

convection  velocity 

shear  velocity 

complex  multiplier  applied  to  the 
element  (Eq.  2.1) 
element  width 

autocorrelation  function  of  the  point  surface 

pressure  under  the  TBL 

Dirac  delta  function 

see  Eq.  A- 10 

channel  gain  deviation 

hydrophone  posit. ion  deviation 

channel  phase  deviation 

t hi 

fractional  amplitude  error  on  the 
e lement 
fluid  density 
reciprocal  wavelength 

rms  channe 1-to-channel  gain  tolerances 
k-space  rms  bandwidth  (Eq  2.7) 

rms  channel -to-channcl  position  tolerances 
rms  channel-to-channel  fractional  amplitude 
rms  channel-to-channel  phase  tolerance 
frequency  ( 2it  times  reciprocal  period) 
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